Aperiodic dynamical decoupling sequences in presence of pulse errors 
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Dynamical decoupling (DD) is a promising tool for preserving the quantum states of qubits. 
However, small imperfections in the control pulses can seriously affect the fidelity of decoupling, 
and qualitatively change the evolution of the controlled system at long times. Using both analytical 
and numerical tools, we theoretically investigate the effect of the pulse errors accumulation for two 
aperiodic DD sequences, the Uhrig's DD (UDD) protocol [G. S. Uhrig, Phys. Rev. Lett. 98, 100504 
(2007)], and the Quadratic DD (QDD) protocol [J. R. West, B. H. Fong and D. A. Lidar, Phys. Rev. 
Lett 104, 130501 (2010)]. We consider the implementation of these sequences using the electron 
spins of phosphorus donors in silicon, where DD sequences are applied to suppress dephasing of the 
donor spins. The dependence of the decoupling fidelity on different initial states of the spins is the 
focus of our study. We investigate in detail the initial drop in the DD fidelity, and its long-term 
saturation. We also demonstrate that by applying the control pulses along different directions, 
the performance of QDD protocols can be noticeably improved, and explain the reason of such an 
improvement. Our results can be useful for future implementations of the aperiodic decoupling 
protocols, and for better understanding of the impact of errors on quantum control of spins. 

PACS numbers: 03.67.Pp, 03.65.Yz, 76.30.-v, 75.10.Jm 
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I. INTRODUCTION 



Mitigating the effect of decoherence is an impor- 
tant problem in the emerging area of quantum in- 
formation processing (QIP)i and other quantum-based 
technologies^Ti^. Dynamical decoupling (DD) is a promis- 
ing tool for this tas k^°'^^ . Originating from the ideas 
that underlie the spin echo effect in nuclear magnetic 
resonance (NMR)^!^^, DD employs a specially designed 
sequence of control pulses applied to the qubits (or the 
central spins) in order to negate the coupling of the cen- 
tral spins to their environment. A variety of DD pro- 
tocols have been developed, analyzed, and employed for 
a few decades in the area of high-resolution NMRMi^. 
It has been suggested later that the decoupling prop- 
erty of these protocols can be used for more general 
purpose of suppressing decoherence and achieving high- 
fidelity quantum operations, thus giving rise to extensive 
theoreticaliiii^"— and experimental22ri£ investigation of 
various novel aspects of DD in the QIP context. Since 
DD requires no feedback and no ancilla resources, and 
is applicable to a wide range of systems, it presents a 
promising tool for lowering the number of quantum er- 
rors beyond the error correction threshold^. 

Performance analysis of the decoupling protocols is of- 
ten based on the Magnus expansion (ME)ii, which is an 
asymptotic cumulant expansion of the evolution opera- 
tor of the system, with the characteristic inter-pulse delay 
(or duration of a single DD cycle, for periodic sequences) 
playing the role of the small expansion parameter. By 
designing the protocols which nullify more and more ex- 
pansion terms, one expects that in favorable experimen- 
tal situations the decoupling fidelity increases. An ex- 
ample of the first-order decoupling protocol is a very ba- 
sic periodic dynamical decoupling (PDD) sequenc o^^'^^ , 
which consists of repeating the basic cycle d — n — d — n, 



where d denotes the free evolution of the system, and tt 
denotes a 180° spin rotation around e.g. x-axis (all anal- 
ysis here and below is performed in the rotating frame, 
i.e. in the coordinate frame rotating around the z-axis 
with the Larmor precession frequencj^) . This sequence 
works for the Ising-type coupling, for instance, when a 
central spin is coupled to the decohering bath via the 
term S^B, where is the z-component of the central 
spin and B is the bath operator. Such terms lead to 
pure dephasing of the central spins in the x-y plane, i.e. 
only transverse components of the central spins are af- 
fected by the bath. In order to eliminate more general 
type of coupling, involving decoherence of all compo- 
nents of the central spin, a two-axis control is needed, 
when the 180° rotation axes alternate between x and y 
directionsii, so the basic period of the PDD sequence 
has a form d — ttx — d — Try — — ttx — d — Try, with 
Trx and Try denoting the Tr rotations around x and y, 
respectively. To increase the decoupling order, eliminat- 
ing the second-order ME terms as well, we need to use 
more advanced sequences than PDD. One possibility is to 
use symmetrized sequences {SDD)^i^^^ having a twice 
longer cycle, consisting of the PDD cycle followed by its 
mirror image. Another approach to achieve higher-order 
decoupling is to use concatenated DD (CDD) protocol, 
which removes the system-bath coupling terms in the ME 
up to the n-th order by nesting the PDD sequence recur- 
sively to itself n times J^ii^iS^ii^ The number of pulses 
in CDD increases exponentially with the concatenation 
level, approximately as 4". Within CDD approach, the 
total decoupling time can be increased either by increas- 
ing the concatenation order, or by periodically repeating 
the fixed-level CDD protocol. 

Another way to improve the decoupling fidelity is 
to optimize positions of the pulses, thus creating a 
family of aperiodic sequences where the inter-pulse de- 
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lays are not necessarily commensurate. In Ref. [l^ 
G. S. Uhrig proposed a sequence (later abbreviated as 
UDD), where optimization of the pulse positions pro- 
vides n-th order decoupling using n pulses^i^. UDD has 
been implemented and tested experimentally in different 
systems^"— i^"—. Theoretical studies have shown that 
UDD is optimal when the noise spectrum of the bath has 
a sharp high-frequency cutof^i^i^, while for the spec- 
tra with soft cutoffs the protocols with periodic struc- 
ture, such as the Carr-Purcell-Meiboom-GillH (CPMG) 
sequence, perform better— i^i^lii^. While UDD employs 
a single-axis control, and therefore suppresses only pure 
dephasing, several UDD-based aperiodic protocols were 
developed to suppress general decoherence by incorporat- 
ing double- axis controls, including the concatenated ver- 
sion of UDD (CUDD)i^, and quadratic DD (QDD) where 
two UDD sequences based on spin rotations about per- 
pendicular axes are nested into each other—. In CUDD, 
the number of pulses also grows exponentially but slower 
than CDD, and in QDD the pulse number required to 
suppress n-order decoherence is of order of 0{n^). 

While DD protocols are usually designed based on the 
assumption of ideal pulses, in reality the accumulation 
of small imperfections of the pulses can severely affect 
the decoupling fidelity and even qualitatively change the 
evolution of the central spins. The influence of the pulse 
errors on the DD protocols has been studied from early 
days of NMRi^iii within the ME (or ME- like) settings, by 
considering the pulse errors as extra terms in the Hamil- 
tonian. However, such an approach is not always satis- 
factory. First, the convergence conditions of ME are not 
always satisfied in the experiments, although the DD may 
still efficiently suppresses decoherence. Second, due to 
its asymptotic nature, the ME can miss the possibility of 
disastrous pulse error accumulation at long times, when 
the number of pulses becomes very large. Very detailed 
studies of the errors introduced by the finite pulse width 
have been performed for CDD^^ and UDDHd*. The ex- 
perimental study of the systematic errors in the rotation 
axis and angle has been done for a number of periodic 
protocols and for UDD^i^, and both error types (finite- 
width and systematic) have been discussed for dynami- 
cally corrected gates^. The theoretical-experimental in- 
vestigation of the long-term accumulation of the system- 
atic pulse errors has been done for several periodic-based 
protocols, PDD, SDD, and CDTM>^. 

In this work, we present a theoretical analysis of the 
aperiodic protocols, UDD and QDD, focusing on the 
long-term accumulation of the systematic pulse errors. 
These errors can be much more devastating than the im- 
perfections associated with the finite width of the pulses. 
We use numerical modeling to analyze the situation when 
the aperiodic protocols with imperfect pulses are used for 
decoupling of the electron spins of the phosphorus donors 
in silicon. This system may present a promising platform 
for scalable solid-state quantum information processing, 
and is a good testbed for studying a number of funda- 
mental issues in quantum dynamics and quantum con- 



trol of solid-state spins^""— . The Si:P system is well 
characterized^"—, DD of the electron spins of P donors 
is efficiently implemented via pulsed ESR technique, and 
the realistic pulse error parameters have been estimated 
earlieriSi^. Comparing the results of our study below 
with the previous works on periodic protocols in the same 
system, we find that the aperiodic sequences provide rea- 
sonable decoupling fidelity, but they are more demand- 
ing to the pulse error magnitude than e.g. CDD. We also 
demonstrate that, by choosing a right set of two-axis con- 
trols, the effect of the pulse errors can be significantly re- 
duced, and explain the reason behind this improvement. 
It is interesting that the right choice for QDD protocol, 
employing rotations around the z and y axes, is a bad 
choice for CDD and PDD, where the rotations around 
X and y should be used^i^. It demonstrates again that 
the analysis of the pulse error accumulation is, in general, 
protocol-specific. 

The rest of the paper is organized as follows. In SecHIl 
we describe the Si:P system and present a model for the 
pulse errors. In Sec. IIIIl we present the analytical and 
numerical results for UDD and QDD. Conclusions are 
given in Sec. |IVl 



II. DESCRIPTION OF THE SYSTEM 

A. The Phosphorus Doped SiHcon System 

The electron spins of P donors in silicon have long re- 
laxation and coherence times^"— , and have been used 
earlier to explore fundamental aspects of quantum spin 
dynamics and DD^i^i^i^i^i^-^. In this study, we 
consider the same experimental situation as described in 
Ref.fi^. The sample is an isotopically purified bulk silicon 
with ^^Si concentration 800 ppm, and very low dop- 
ing density of phosphorus, so that the coupling between 
different P donor spins can be neglected at the relevant 
timescales. A large static (quantizing) magnetic field is 
applied along the z-axis to the system, and the tempera- 
ture is sufficiently low (8 K in the experiments reported 
in Ref. [ioi ) . The relaxation time Ti of the P electron spin 
is very long, so the longitudinal relaxation is neglected 
here. Similarly, the coherence time T2, defined as the 
spin echo decay time, is long (several milliseconds, lim- 
ited by instantaneous diffusion) and will also be neglected 

Thus, we consider only pure dephasing of the P 
electron spins. One important dephasing channel is 
the Ising-type hyperfine coupling to the ^^Si nuclear 
spins^"— i^i^. Due to rather large localization radius 
of the P donor electron in silicon, the P electron spin 
couples to a large number (bath) of ^^Si nuclear spins. 
For an isotopically purified sample considered here, the 
interactions between different ^^Si nuclear spins are small 
and can be neglected, so we are dealing with the static 
nuclear spin bath. Moreover, since the P donors are well 
separated, the ^^Si nuclei which are appreciably coupled 
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to one P electron spin interact very weakly with other 
donor centers, so each donor can be considered as a sep- 
arate central spin with its own static spin bath. An- 
other important contribution to dephasing is the inho- 
mogcneous broadening due to quasistatic fluctuations of 
the quantizing magnetic field. Therefore, in the frame 
rotating with the frequency at the center of the ESR line 
of the P electrons, the Hamiltonian for a P electron spin 
coupled to the bath is 

n = (1) 



B. The analytical model of pulse errors 

For the experimental situation under consideration the 
systematic errors in the pulse parameters (rotation angle 
and the direction of the rotation axis) are much larger 
than the errors associated with the finite pulse width. 
Therefore, we treat each 7r-pulse as an instantaneous but 
imperfect spin rotation. For a nominal 7r-pulse about the 
X axis (the ttx pulse) the unitary operator describing the 
spin evolution is 

C/x = expH(^ + e:.)(S-n)] (3) 



where S*^ is the z-component of the electron spin opera- 
tor, and B is the total random static offset field created 
by the bath of ^^Si nuclei and the fluctuations of the 
external fieldi^i^. Experiments show that the ESR line- 
shape is Gaussian with the width of 50 mG, so that the 
distribution function of B is 



where is the error in the rotation angle, and n 



P{B) 



V2^ 



exp 



B2 



(2) 



with 6 = 50 mG. Note however, that our results are not 
very sensitive to the specific shape of the distribution. 
The hyperfine coupling of the P electron to the donor's 
own '^^P nuclear spin (/ = 1/2) is larg o^^'^^ , of order of 
100 MHz, so that the two hyperfine lines are well sepa- 
rated, and we assume that only one line is excited. Then 
the on-site hyperfine coupling to '^-'^P nucleus just shifts 
the ESR frequency, and we neglect this interaction. 

As is standard in NMR and ESR, the experimental 
temperature is much larger than the Zeeman energy of 
the electron spins in the quantizing field. Thus, we de- 
scribe the state of a P spin within the high-temperature 
approximation, as is customarily done in NMR/ESR 
theorjii^, using the density matrix p cx 1 -t- aS'^ , where 
a is small. Since the Ti time is extremely large, and 
the identity matrix is not affected by the unital evolu- 
tion, we can neglect it, and consider the electron spin as 
being in a pure state with = 1/2 (so-called pseudo- 
pure state). Other initial (pseudo)pure states, e.g. along 
the X and y axes, can be prepared by the 7r/2 pulses ap- 
plied along the y and —x axes, respectively. For pure 
initial states performance of the decoupling can be con- 
veniently quantified by the input-output fidelity, i.e. by 
the overlap between the initial and the final (reduced) 
density matrices of the system. By virtue of the relation 
p{t) = ■| + (S'^)crx + (5'^)cry + (S'^)o-2, the average spin pro- 
jections (5^), (5*^), and (S*^) can be equivalently used as 
measures of the decoupling performance. Note that the 
angular brackets here denote both quantum-mechanical 
averaging (including trace over the bath states) and the 
averaging over the ensemble of different P electron spins. 
Below, we use the rescaled fidelity, which for the initial 
state along the axis a = x,y^z is defined as Fa — 2{Sa)- 



1 — riy — n1,ny,nz) is the unit vector along the ro- 
tation axis, where the small parameters Uy and char- 
acterize the deviation of the actual rotation axis from the 
X axis. Similarly, the rotation operator for a nominal ny 
pulse is 

C/y = exp[-i(7r-fey)(S-m)] , (4) 

where ey is the rotation angle error, and lii = 
(toj;, -^1 — — , niz) is the actual rotation axis, with 
small parameters nix and niz ■ During the delay between 
neighboring pulses, the evolution operator is simply 



Ud{T) = cxp[~iBS'T] 



(5) 



where r is the duration of the delay between the pulses 
(e.g., for periodic sequences studied in Ref. equal to 
11 fis). 

To model the systematic errors, we take into account 
that the resonant ac field Bp, which rotates the spins 
during pulses, is not homogeneous over the sample. We 
assume for simplicity that this field varies only along 
one spatial axis, denoted as I, and the sample is located 
between I = +d and I = —d. Wc also assume that 
the sample is placed optimally with respect to the field, 
with the sample center I — located at the maximum 
of Bp. Then the spatial dependence of Bp has a form 
Bp{l) = Bp + ABp[l — 3P/d^], where Bp is the average 
magnitude of the ac field over the sample and ABp char- 
acterizes the ac field inhomogeneity. If the pulse width 
tp is adjusted to give Bptp = tt, then the spins located 
at different parts of sample will undergo rotation by dif- 
ferent angles. For ttx pulse, the resulting rotation angle 
errors ex will have distribution 



P(e,) = (l/2eo)[3(l-e./eo)] 



-1/2 



(6) 



where — 2eo < e < coi and cq characterizes the magnitude 
of the rotation angle errors. In a similar way, we assume 
that the z-component of the rotation axis Uz has the same 
distribution, and its magnitude is characterized by the 
parameter ng. Wc also assume that the pulse errors for 
TTx and tty pulses are the same, taking — and = 
Uz ■ The pulse error parameters rrix and ny have different 
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nature: they characterize the phase of the pulse, and can 
be adjusted very precisely, so we assume them to be zero, 
following Refs. EolFll . Thus, within this model, all pulse 
errors are characterized by only two parameters, eg = 0.3 
and Tin = —0.12. These values have been determined in 
Refs. Eollsil . and have reasonable magnitude. In spite 
of utter simplicity, this model reproduces (or, at least, 
mimics) well most features observed in experiments. 



III. NUMERICAL AND ANALYTICAL 
RESULTS 

A. Performance of UDD 

Uhrig's DD protocol, UDD, is based on a single-axis 
control, i.e., all spin rotations arc performed about the 
same nominal axis, which we take as x. For UDD of 
level denoted as UDD-£, the iix pulses arc applied at 
timesi^ 



t sin 



2£ + 2 



(7) 



where t is the total evolution time, and j = 1, 2, . . . , ^ for 
even (, and j = 1, 2, . . . , ^ + 1 for odd The simulation 
results show that UDD of even levels and UDD of odd 
levels behave differently. Fig. [1] shows the decoupling 
fidelities for UDD-£ as functions of total evolution time 
for £ = 2, 20, and Fig. [2] shows the resuhs for UDD-€ 
with ^ = 3 and 19. 
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FIG. 1: (Color online) . Decoupling fidelities of UDD-2 (a) and 
UDD-20 (b) as functions of the total evolution time for initial 
states Sx (red empty circles, red lines), Sy (green triangles, 
green lines), and Sz (blue soUd circles, blue lines). 



We consider fidelities for three different initial states of 
the central spins, directed along x, y, and z axes. They 
characterize how well the corresponding spin components 
are preserved. For UDD-2, which is the same as a single 
cycle of CPMG, the spin component Sx is preserved with 
fidelity close to 1, while the fidelities for the other two 
components are affected by the pulse errors. For long evo- 
lution time, Fy and saturate at values ~ 0.89. Strong 
dependence of the fidelity on the initial state has been no- 
ticed and analyzed before for periodic protocols, and it 
is not unexpected that aperiodic sequences demonstrate 



this feature as wel^^^'^^. For I = 20, Sx is well pre- 
served, demonstrating only a slow decay, while Fy and Fz 
have initial values close to zero, and saturate to a value 
around 0.5. Such low values even for t = are caused 
purely by a catastrophic accumulation of errors of 20 ttx 
pulses. 

Let us consider the level-2 UDD as an example in order 
to examine the long time behavior. The evolution oper- 
ator, up to second order in the pulse errors, and using 
Uy = 0. is 



(8) 



■ cos 4^ + 2n, sin 4^ and 9-, = e^nz cos 



with 9.j; — ^v^o ^ I ^,„z^iii ^ LOii^ — ^x'^'Z^^'^ 2 

{n1 — sin Although the qualitative features of the 
protocol can be analyzed using only the first-order ex- 
pression for quantitative analysis of the fidelity re- 
quires the second-order terms in this evolution operator. 
Taking into account that e^, Uz, and B are independent 
variables, and that (e^:) = (n^) — 0, the expression for 
Fy becomes 



F„ = l-20i 



(9) 



l-2[(4)(cos^f >+4(n^)(sin2^)]. 



Noticing from Eq. 
O.Shq, wc arrive at 



that 



O.Seq, and 



lim Fy = l- 0.8(4 + 4?io) = 0- 



(10) 



This is close to the value 0.89 which we see in the sim- 
ulations shown in Fig. [l]for UDD-2. The analysis above 
relies on the fact that B is time-independent, so in ex- 
periments such saturations will be replaced by decays at 
different rates, depending on the bath dynamics. There- 
fore, although our results extend until t — 2000 /is, at 
these times the instantaneous diffusion and the internal 
bath dynamics should be taken into account. We neglect 
these effects here since they are beyond the scope of our 
present study, and require a separate focused research 
effort. 

For UDD of higher level, similar behavior is observed. 
In Fig. [TJ the fidelities for UDD-20 are generally smaller 
than for UDD-2, which reflects more serious accumula- 
tion of the pulse errors. 

An interesting feature, seen in Fig. [U is that at longer 
times, when dephasing becomes noticeable, the fidelities 
may increase, as if decoherence counteracts the pulse er- 
rors. To understand this, let us consider UDD-2 in the 
situation when the errors arc absent, and only the 
rotation angle error ex is nonzero. If a spin is prepared 
along the y axis, and is subjected to a sequence of imper- 
fect TTx pulses, then its y-component rotates around the 
X-axis, and accumulates the rotation angle errors. But if 
the random field B moved the spin away from the y-axis 
towards the x-axis by the angle x = Bt/A, then the spin 
component along the x axis is not affected by the rotation 
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angle errors, and the spin is refocused with better pre- 
cision. Correspondingly, the fidelity for a fixed value of 
B demonstrates a modulation proportional to — cos^ x, 
with maximum fidelity achieved at x = 7r/2. 




t (us) t (us) 



FIG. 2: (Color online) . Decoupling fidelities of UDD-3 (a) and 
UDD-19 (b) as functions of the total evolution time for initial 
states Sx (red empty circles, red lines), Sy (green triangles, 
green lines), and Sz (blue solid circles, blue lines). 

For UDD of odd level i = 2n — 1, the number of pulses 
is 2n since there is a tt pulse at the end of the whole 
evolution. As shown in Fig. [51 for a given UDD-£, when 
the total evolution time is short, Sx is better preserved 
than the other two spin components. At long time, 
and Fy saturate at similar values, while F^ saturates at 
a lower value. The process of saturation is much slower 
for £=19 than for £ = 3. 

To gain insight into the behavior of odd-level UDD, 
we consider UDD-3. The evolution operator for this se- 
quence, up to first order in the pulse errors, is 

J/g™ = 1 - iOnCTx - i'HnCy (H) 

with 

0n = y[l-H2cOs(BTi)+COs(B(T2-Ti))] 

+ n^[2sin(BTi) + sin(B(T2 -ri))] (12) 

and 

Vn = - 2cOs(Bti) + COS (B(t2 - Ti))] 

- |-[2sin(STi) - sin {B{t2 - ti))] , (13) 

where ri and T2 refer to the durations of the first and 
second inter-pulse delays. In general, for quantitative 
calculation of the fidelities, we need also to consider 
the second-order terms, but it can be shown that for 
UDD-3 their contribution is zero. At long times, when 
(sin(i3Ti)) = (sin (i3(r2 — n))) ~ (and similarly for 
cosines), the fidelity Fy acquires a form 

Fy = l-2{0l) 

= l-2{hnl+'l + Cl^nl)(Acos\BT,) 

+ cos" {B[t2-t^)))) (14) 



and wc finally obtain 

lim Fy = l- b{nl) - = 0.82, (15) 

which coincides with the result of the numerical simula- 
tions. 

We now examine the fidelities at i — > 0. Since we treat 
pulses as instantaneous rotations, this situation should 
be understood as the limit where the total evolution time 
is much smaller than the dephasing time T2 = 1/b, but 
much larger than the total width 2ntp of all pulses. Decay 
of fidelity in this case arises only from accumulation of 
the pulse errors. The evolution operators for UDD of 
both even {£ = 2n) and odd {£ = 2n — 1) levels have the 
same expression 

limC/r = i-ini - me^a,) , (16) 
t-yO 

corresponding to a spin rotation about the x axis with 
the rotation angle proportional to n. Small differences in 
the rotation angles for different spins are scaled by n. For 
large n the spins become almost uniformly distributed in 
the y-z plane, hence Fy and F^ are close to zero, while 
Fx is close to 1, as seen in Figs. [T] and [2] for £ = 20 and 
£= 19. 



B. Performance of QDD 

In order to suppress general decoherence, DD sequence 
has to use two-axis control. An extension of UDD, 
the quadratic DD (QDD) has been suggested for this 
purpose^^. QDD is composed of two nested UDD se- 
quences, with inner and outer sequences employing the 
rotations about two diff'erent perpendicular axes. A QDD 
sequence of order £, which we denote as QDD-^, can sup- 
press decoherence to £-th order terms with {£+ 1)^ intcr- 
pulse intervals. 



1. QDD based on ttx and Try pulses 

We first consider the QDD sequence with the ttx pulses 
in the outer hierarchical level, and Try pulses in the inner 
hierarchical level, examining the effects of error accumu- 
lation. 

Like UDD, QDD of even and odd levels also behave 
differently. For £ = 2?i — 1, the sequence of QDD-£ is 

UDDf VJ-t^x-UDD^ V2)-7rx • • • -UDD^ Vf+i)-^x 

(17) 

where E^tl = t and UDD^'^Vi) denotes UDD-£ 
based on Try pulses with evolution time Tj . The division 
of the total time t into intervals Tj follows the same rule 
as UDD, i.e., Tj = tj—tj_i, with tj given by Eq. ([7]). The 
number of pulses in QDD with £ = 2n-l is {£+l){£ + 2). 
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For £ = 2n the sequence QDD-f is 

UDDP(ti)-^x- • • • -UDD(^'(T,)-7rx-UDD(^)(r,+i) 

(18) 

and the number of pulses is £{£ + 2). 
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FIG. 3: (Color online). Decoupling fidelities as functions of 
the total evolution time for initial states Sx (red empty cir- 
cles, red lines), Sy (green triangles, green lines), and Sz (blue 
solid circles, blue lines), (a) QDD-3 (20 pulses), (b) QDD-4 
(24 pulses). 

Considering the limit i — )■ 0, we have for even and odd 
levels 

limC/2T = I - in{e^;(Tx + f^yO-y) (19) 
limC/r-i = (-!)"(! -*ne,a,) (20) 

That is, for short total evolution time, in a sequence 
QDD-(2n), the dominant pulse error is e^, while in QDD- 
(2n — 1) both and Cy contribute in first order. The 
difference is due to an extra ttx pulse at the end of the 
odd-order QDD. 

Simulation results for QDD-3 are shown in Fig. [3] (a) . 
Fidelities for all three spin components exhibit long-time 
saturation. The saturation value for Sx is higher than 
for the other two components. However, Fx in QDD-3 
exhibits initial decay (even for t = 0), instead of being 
preserved as expected from Eq. ([20)1 . This is because for 
the system under consideration, the error parameter eg = 
0.3 is pretty large, and the pulse errors may appreciably 
contribute in higher orders. Taking QDD of = 3 as an 
example, in the limit < — > 0, the sequence is a fourfold 
repetition of the unit 

(7rY-7rY-7i'Y~7rY)~'''x • (21) 

Expanding the corresponding evolution operator up to 
second order in the pulse errors, and taking into account 
iTT-x = TLy = 0, we find that 

U = {l- 2el) - 2iex + 2ie^(2ej, -f- n,) . (22) 

Thus, the rotation axis deviates from the x axis towards 
the z axis, and the spin state along the x axis is not ex- 
actly parallel to the rotation axis, resulting in the initial 
drop of Fx- 

Fig- El (b) shows the simulation results for QDD-4, with 
24 pulses. The saturation values in a descending order 



are Fy, F^ and F.. The initial values for F^. and Fy are 
the same, which can be explained by Eq. taking into 
consideration that Cx = ey in our model. 

As we see, the QDD sequences are rather sensitive to 
the pulse errors. While the saturation at long times is 
seen for all spin components, the initial fidelity loss is 
noticeable. 



2. QDD based on nz and ny pulses 

In order to further understand the effect of pulse er- 
rors, we briefly consider another QDD sequence with wz 
pulses in the outer hierarchical level and Try pulses in the 
inner one. We denote such a sequence as QDD(ZY). To 
produce nz pulses in ESR experiments, the pulses ttx 
and TTy are applied back-to-back^, and our simulations 
reproduce this situation, taking into account the ttx and 
TTy pulse errors. Fig. [?] shows the simulation results for 
QDD(ZY) of levels 3 and 5 (two upper panels), and 2 
and 4 (two lower panels), with 20 and 42, and 8 and 
24 pulses, respectively. While QDD(ZY) of even levels 
provide no obvious improvement compared to QDD se- 
quences based on ttx and Try pulses, the QDD(ZY) pro- 
tocols of odd level preserve all three spin components, 
with high saturation values at long times. 

To understand this effect, let us inspect the structure 
of the QDD(ZY) sequence. As seen from Eq. (|17p . since 
UDD of odd level has a Tr-pulse at the end of the se- 
quence, after nesting this Try pulse is followed by a tt^ 
pulse without delay, and these two adjacent pulses are 
equivalent to a (imperfect) ttx pulse. The sequence of 
QDD(ZY) thus has a structure 

{Tl,l-''^Y-Tl,2-TrY-Tl,3-'!TY-Tl^4)-1Tx- 
(T2,l-TrY-T2^2-7rY-T'2,3-7rY-T2,4)-7rx- 
('''3,l-TrY-T3,2-7rY-T'3,3-7rY-T'3,4)-7rx- 

(T4,i-7rY-T'4,2-7rY-T'4,3-7rY-T'4,4)-7rx (23) 

where Tij denotes the j-th delay in the i-th UDD in the 
whole sequence. Therefore, QDD(ZY) of odd level is ac- 
tually based on pulses about x and y axes, and has a 
structure similar to the periodic sequences based on x 
and y rotations, which are known to be very robust with 
respect to the pulse error o'^'''^^i"'^°'^^i^^'^"'^ . Such a struc- 
ture can not be achieved with QDD based on ttx and Try 
pulses. At even level, the x,?/-based QDD lacks a final 
TTx pulse, see Eq. [13 At odd level, the x, y-based QDD 
has even number of Try pulses sandwiched by ttx pulses, 
see Eqs. [T7] and [HB while for robustness we need odd 
number of y-pulses. 

In the limit t — )■ 0, to first order in the pulse errors, 
the evolution operator for QDD(ZY) of level = 2n — 1 
is 

Utln-i = (-1)"[1 + 2m(m, + ny)a,]. (24) 

Since for our system mx ^ ny = 0, the DD performance 
is determined by accumulation of the higher-order errors. 
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Therefore, for QDD(ZY) of odd level, the dominant 
pulse errors are the in-plane components of rotation axis 
errors. For the system under consideration such errors 
are negligible, and QDD(ZY) is a good choice for deco- 
herence suppression for all quantum states. 
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QDD, using the phosphorus electron spins in silicon as 
an example of experimental implementation. We show 
that the decoupling fidelity strongly depends on the ini- 
tial state. At long times, we observe saturation of fideli- 
ties. UDD and QDD of even and odd levels are found to 
perform differently. 

UDD protocol is based on a single-axis control, so 
the decoupling fidelity is highest for the spin component 
along the control axis. QDD preserves the three spin 
states with closer fidelities. In particular, QDD of odd 
level based on -kz and Try pulses is very robust with re- 
spect to the pulse errors, due to its special structure. Our 
results can be useful for experimental implementations of 
the aperiodic decoupling protocols, and for deeper under- 
standing of the influence of errors on quantum control of 
spin systems. 



FIG. 4: (Color online). Decoupling fidelities for QDD(ZY) as 
functions of the total evolution time for initial states Sx (red 
empty circles, red lines), Sy (green triangles, green lines), and 

Sz (blue solid circles, blue lines), (a) level 3 QDD(ZY), (b) Acknowledgments 
level 5 QDD(ZY), (c) level 2 QDD(ZY), (d) level 4 QDD(ZY). 



IV. CONCLUSION 

We have studied accumulation of the pulse errors for 
two aperiodic dynamical decoupling sequences, UDD and 



We would like to thank K. Khodjasteh and L. Viola for 
useful discussions, and A. M. Tyryshkin and S. A. Lyon 
for useful discussions and stimulating suggestions. Work 
at the Ames Laboratory was supported by the Depart- 
ment of Energy - Basic Energy Sciences under Contract 
No. DE-AC02-07CHII358. 



^ M. A. Nielsen, I. L. Chuang, Quantum Computations 
and Quantum Information (Cambridge Univ. Press, Cam- 
bridge, 2002). 

^ L. Childress, J. M. Taylor, A. S. Sorensen, and M. D. 

Lukin, Phys. Rev. Lett. 96, 070504 (2006). 
3 J. A. Jones et al, Science 324, 1166 (2009). 
* P. Cappellaro, J. Emerson, N. Boulant, C. Ramanathan, 

S. Lloyd, and D. G. Cory, Phys. Rev. Lett. 94, 020502 

(2005). 

^ J. M. Geremia, J. K. Stockton, and H. Mabuchi, Phys. 

Rev. Lett. 94, 203002 (2005). 
® J. M. Taylor, P. Cappellaro, L. Childress et al., Nature 

Phys. 4, 810 (2008). 

C. Degen, Appl. Phys. Lett. 92, 243111 (2008). 
^ G. Balasubramanian, I.-Y. Chan, R. Kolesov et al., Nature 

455, 648 (2008). 
^ G. de Lange, D. Riste, V. V. Dobrovitski, R. Hanson, 

arXiv: 1008.4395. 
° L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
^ L. Viola, S. Lloyd, and E. Knill, Phys. Rev. Lett. 83, 4888 

(1999). 

^ E. Hahn, Phys. Rev. 80, 580 (1950). 

^ G. P. Slichter, Principles of Magnetic Resonance (Springer, 
Berlin, New York, 1990). 



U. Haeberlen, High Resolution NMR in solids: Selective 

Averaging (Academic Press, New York, 1976). 

B. G. Gerstein and C. R. Dybowski, Transient Techniques 

in NMR of Solids (Academic Press, Orlando, 1985). 

K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 

180501 (2005). 

Ren-Bao Liu, Wang Yao, L. J. Sham, New J. Phys. 9, 226 
(2007). 

W. X. Zhang, V. V. Dobrovitski, L. F. Santos, L. Viola, 
and B. N. Harmon, Phys. Rev. B 75, 201302(R) (2007). 
G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 
^° L. Cywinski, R. M. Lutchyn, CP. Nave, S. Das Sarma, 
Phys. Rev. B 77, 174509 (2008). 

L. Cywinski, W. M. Witzel, and S. Das Sarma, Phys. Rev. 
B 79, 245314 (2009). 
2^ G. S. Uhrig, New J. Phys. 10, 083024 (2008). 

W. Zhang, N. P. Konstantinidis, V. V. Dobrovitski, B. N. 
Harmon, L. F. Santos, and L. Viola, Phys. Rev. B 77, 
125336 (2008) 

W. Zhang, N. P. Konstantinidis, K. ALHassanieh, V. V. 
Dobrovitski, J. Phys.: Gond. Matter 19, 083202 (2007). 
^■■^ J. R. West, B. H. Fong and D. A. Lidar, Phys. Rev. Lett 
104, 130501 (2010). 

K. Khodjasteh, D. A. Lidar, and L. Viola, Phys. Rev. Lett 



8 



104, 090501 (2010). 

S. Pasini and G. S. Uhrig, Phys. Rev. A 81, 012309 (2010). 

E. Fraval, M. J. Sellars, and J. J. Longdell, Phys. Rev. 

Lett. 95, 030506 (2005). 
23 J. J. L. Morton et al, Nature (London) 455, 1085 (2008). 
^° D. Press, T. D. Ladd, B. Zhang, Y. Yamamoto, Nature 

456, 218 (2008). 

H. Uys, M. J. Biercuk, and J. J. Bollinger, Phys. Rev. Lett. 
103, 040501 (2009). 
^2 M. J. Biercuk, H. Uys, A. P. Vandevender, N. Shiga,W. 
M. Itano, and J. J. Bollinger, Nature (London) 458, 996 
(2009). 

J. Du, X. Rong, N. Zhao, Y. Wang, J. Yang and R. B. Liu, 
Nature 461, 1265 (2009). 

V. V. Dobrovitski, G. de Lange, D. Riste, R. Hanson, Phys. 
Rev. Lett. 105, 077601 (2010) . 
3^ B. Naydenov et al, arXiv:1008.1953v2 (2010). 

G. A. Alvarez, A. Ajoy, X. Peng, and D. Suter, Phys. Rev. 
A 82, 042306 (2010). 

G. de Lange, Z. H. Wang, D. Riste, V. V. Dobrovitski, R. 
Hanson, Science 330, 60 (2010). 

C. A. Ryan, J. S. Hodges and D.G. Cory, Phys. Rev. Lett. 

105, 200402 (2010). 

^3 H. Bluhm et al, ArXiv:1005.2995vl (2010). 

A. M. Tyryshkin, Zhi-Hui Wang, Wenxian Zhang, E. 
E. Haller, J. W. Ager, V. V. Dobrovitski, S. A. Lyon, 
arXiv:1011.1903v2. 

"1 L. F. Santos and L. Viola, Phys. Rev. A 72 062303 (2005). 
"2 K. Khodjasteh and D. A. Lidar, Phys. Rev. A 75, 062310 

(2007) . 

B. Lee, W. M. Witzel, and S. Das Sarma, Phys. Rev. Lett. 
100, 160505 (2008). 

W. Yang and R.-B. Liu, Phys. Rev. Lett. 101, 180403 

(2008) . 

Ashok Ajoy, Gonzalo A. Alvarez, Dieter Suter, 



arXiv:1011.6243v2. 
■"^ G. S. Uhrig, Phys. Rev. Lett. 102, 120502 (2009). 

R. W. Vaughan, D. D. EUenian, L. M. Stacey, W.-K. Rhim, 

J. W. Lee, Rev. Sci. Instr. 43, 1356 (1972). 

G. S. Uhrig, S. Pasini, New J. Phys. 12, 045001 (2010). 
^3 M. J. Biercuk, H. Uys, A. P. VanDevender, N. Shiga, W. 

M. Itano, and J. J. Bollinger, Phys. Rev. A 79, 062324 

(2009). 

K. Khodjasteh, L. Viola, Phys. Rev. A 80, 032314 (2009). 
'^^ Zhi-Hui Wang, Wenxian Zhang, A. M. Tyryshkin, S. A. 

Lyon, J. W. Ager, E. E. Haller, and V. V. Dobrovitskil, 

arXiv:1011.6417vl. 
^2 B. E. Kane, Nature 393, 133 (1998). 

■''^ A. M. Tyryshkin et al, J. Phys.: Cond. Matter 18, 783 
(2006). 

^'^ A. R. Stegner et al, Nature Phys. 2, 835 (2006). 

A. Morello et al, arxiv: 1003.2679 (2010). 
'^^ A. Honig and J. Combrisson, Phys. Rev. 102, 917 (1956). 

G. Feher and E. A. Gere, Phys. Rev. 114, 1245 (1959). 
''^ A. M. Tyryshkin, S. A. Lyon, A. V. Astashkin and A. M. 

Raitsimring, Phys. Rev. B 68, 193207 (2003). 
'^^ R. C. Fletcher, W. A. Yager, G. L. Pearson, F. R. Merritt, 

Phys. Rev. 95, 844 (1954). 
^° R. de Sousa and S. Das Sarma, Phys. Rev. B 68, 115322 

(2003). 

W. M. Witzel, R. de Sousa, and S. Das Sarma, Phys. Rev. 
B 72, 161306(R) (2005). 

W. M. Witzel, M. S. Carroll, A. Morello, L. Cywiriski, S. 
Das Sarma, Phys. Rev. Lett. 105, 187602 (2010). 
^''^ A. Abragam, Principles Of Nuclear Magnetism Clarendon, 
Oxford, England, 1961. 

T. GuUion, D. B. Baker, M. S. Conradi, J. Mag. Res. 89, 
479 (1990). 



